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Abstract 

A matricial Darboux operator intertwining two one-dimensional 
stationary Dirac Hamiltonians is constructed. This operator is such 
that the potential of the second Dirac Hamiltonian as well as the cor- 
responding eigenfunctions are determined through the knowledge of 
only two eigenfunctions of the first Dirac Hamiltonian. Moreover this 
operator together with its adjoint and the two Hamiltonians generate a 
quadratic deformation of the superalgebra subtending the usual super- 
symmetric quantum mechanics. Our developments are illustrated on 
the free particle case and the generalized Coulomb interaction. In the 
latter relativistic counterpart of shape-invariance is observed. 
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1 Introduction 



In quantum mechanics, the Schrodinger equations which can be solved by 
analytic methods exclusively are rather exceptional. Therefore the methods 
being able to enlarge the number of such equations have attracted much at- 
tention in recent as well as less recent literature. Three of them still remain 
very popular: the Darboux transformations Q elaborated in 1882 within 
the mathematical framework of Sturm-Liouville differential equations, the 
factorization method introduced by Schrodinger H] in 1940 and more re- 
cently the so-called super symmetric quantum mechanics Q. All of them are 
more or less based on the same following ideas. 

Let us consider the following Schrodinger Hamiltonian 

Hq = --^ + Vq{x), xe^oixeW^ (1) 
which can be factorized as follows 

+ a, a = constant (2) 

with 

Then the eigenfunctions of the isospectral (up to the eventual creation or 
loss of one energy) Hamiltonian Hi defined by 

d? 

Hi = W + a = -—^ + Vi{x) (4) 
ax^ 

are obtained through the application of L to the eigenfunctions of Hq as it 
is clear from the so-called intertwining relation 

LHq = HiL. (5) 

Thus a new Schrodinger Hamiltonian Hi has been constructed and it is 
exactly solvable if Hq is. 

We remark that the Darboux transformation has two particular features 
compared to |^, First, the potential W{x), written as 

W{x) = _lMoM, (6) 
dx 

can be constructed from a bounded or unbounded non-vanishing eigenfunc- 
tion ipo{x) of Hq (with eigenvalue a). Second, the Darboux operator L can 
be extended to higher order ||l], |^ . 

Here we shall ask for the same kind of developments in the relativis- 
tic context that is to say search for the operator L intertwining two one- 
dimensional Dirac Hamiltonians. A partial answer has already been given in 
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and through the super symmetrical features of specific Dirac Hamilto- 
nians. Another one can also be found in 0, |8| where a relativistic Darboux 
transformation has been considered but for pseudoscalar potentials only. In 
the following we will not limit ourselves to such a context and will give, in 
Section the extended intertwining operator L corresponding to a general 
self-adjoint potential. This operator is constructed from two (known) so- 
lutions of the initial Dirac equation and gives rise to new exactly solvable 
Dirac equations. Moreover, in Section ^, we will convince ourselves from 
this operator L and its adjoint that the underlying superstructure in the 
relativistic context is a quadratic deformation of the sqm(2) superalgebra, 
the latter being, as well known [^, the subtending superalgebra of the (non- 
relativistic) supersymmetric quantum mechanics. Finally, in Section ^ we 
will illustrate our statements on two examples: The free particle case and 
the generalized Coulomb interaction. For the latter, we observe the rela- 
tivistic counterpart of the so-called shape-invariance [^, i.e., only the values 
of the parameters introduced in the expression of the potentials change. 



2 Intertwining operator for the Dirac equation 

Let us start with the following one-dimensional Dirac Hamiltonian 

ho = ia2-, — 1-^0(3^), X G M or X € Mrt (7) 
ax 

where (T2 is the usual two-by-two Pauli matrix and vq is real and symmetric, 



v,{x) = [ T r • (8) 

V ^12(2;) ^^22 (2;) J 

We assume here that /iq is a known exactly solvable Hamiltonian; in other 
words, all its eigenfunctions, the two-component spinors ip{x), as well as the 
corresponding energies are analytically determined. Let us now search for a 
matricial operator L satisfying the intertwining relation similar to @, i.e., 

Lho = hiL (9) 

with ^ 

/ii = «o-2— + (10) 

vi{x) being at this level the unknown real and symmetric potential. The 
simplest operator L we can consider is 
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where A and B are two-by-two matrices with x-dependent entries. The 
relations @ and (11) give the following system 



[A, CT2] = 0, 

[B, (T2] — iAvQ + iviA — a2Ax = 0, 
Avqx + Bvo - viB - ia2B^ = 0, 
the notation Ax meaning here 



(12) 

(13) 
(14) 



dA 

dx 



( dAr^ 
dx 

dA^i 
\ dx 



dAu \ 

dx 

dA22 

dx / 



The condition (12) is equivalent to ask for An = A22 and A12 = —^21- The 
constraint (O) enables us to fix the potential difference Av = vi — vq 



Av = {Avo - vqA + i[B, CJ2] - i(J2Ax)A ^ 



(15) 



up to the assumption of the existence of A'^. Finally from ( |T4[ ) we can 
obtain the matrix or in a simpler way a defined through B = An. Indeed, 
Eq. (|T4|) then reads 



(vq - i(J2(T)x + [(T, vq] +i[a2, a]a = 0. 



(16) 



We recognize a matrix analogue of the Riccati equation. It can be linearized 
through the substitution 

a = -UxU~^ (17) 

in order to become 

[u~^{vou + ia2Ux)]x = (18) 
which after integration leads to 

hou = ia2Ux + vqu = uX, (19) 

the matrix A being the constant of integration. 

This equation (19) is thus formally speaking an ordinary Dirac one up to 
the fact that the solution u is not a spinor anymore but a two-by-two matrix 
while the usual energy E has also been replaced by a two-by-two matrix A. 

The next step is to find a convenient u that is a solution of (jl^) being real 
(and invertible) in order to ensure the self-adjointness of vi through ([T5|). It 
is ensured in a straightforward manner if 



(ni,U2) , A = diag(ei,e2) 



(20) 
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with the spinors ui and U2 being eigenfunctions (not necessarily bounded) 
of the Dirac Hamiltonian ho 



L = A{ — 1 (22) 



hoUj = EjUj, j = 1,2. (21) 
Having found u, the operator L given in Eq. (11) or 

d 

dx 

as well as the new potential vi (see Eq. (p^)) 

vi= A {vq + i[cr, 0-2] - ia2A-^A^) A'^ (23) 

are now fixed up to the determination of A. This matrix keeps arbitrari- 
ness: all one knows is that it has to commute with 02- For simplicity and 
comparison with the non-relativistic context, we put A equal to the identity 
matrix. Eqs. ( p2|) and ( p3|) are then simplified as follows 

L = ^-u^u-\ (24) 

vi=Vo + i[cr,a2]. (25) 

These results are the relativistic analogues of the usual Darboux transfor- 
mation. We now give another expression in what concerns vi, particularly 
useful for applications. Indeed from (p^ we have 



and therefore 



I.e., 



f = —UxU ^ = i(T2U\u ^ — i(T2V0 (26) 

vi = a2Voa2 + uXu^^ — a2uXu^^a2, (27) 



£i - £2 f di d2 \ 

where di = uuU22 + 1^12^^211 ^2 = ^^21^^22 — wiiiii2, with Uij corresponding 
to the element of the matrix u being at the crossing of the i*'^ lign and the 
j^^ column. 

Let us close this Section by noticing that, by definition, the operator L 
has a non-trivial kernel since kerL = u. This implies that the action of L to 
an eigenspinor of ho corresponding to an eigenvalue different from ei and 82 
will give rise to an eigenspinor of hi. The eigenspinors of hi related to the 
eigenvalues ei and £2 will be obtained through v = {u^)^'^, that is hiv = vX. 
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3 Factorization properties of Dirac Hamiltonians 
and second order supersymmetry 

Let us here consider in addition to L given in Eq. (p4|), its adjoint defined 
by 

L^ = -^- iuxu'')l (29) 

It satisfies an intertwining relation similar to Eq. 

L^hi = hoLl (30) 

This relation means that the operator realizes the transformation in 
the opposite direction, i.e., the application of to the eigenspinors of hi 
gives us the eigenspinors of /iq- The operator L^L is thus such that applied 
to the eigenspinors of ho, it gives back these eigenspinors. By definition, 
this is nothing but the fact that L^L is a symmetry operator of the initial 
Dirac equation hQip = Eip. Since we limit ourselves to the one-dimensional 
stationary context, this implies that L is a function of ho. Moreover 
because is a second order differential (matricial) operator while ho is 
of the first order, L^L is in fact a polynomial of second order in ho. More 
precisely, after tedious calculations, one can be convinced that 

L^L = {ho-ei)iho-e2) (31) 
while a similar result holds for LL^ 

W = {hi-ei){hi-e2). (32) 
If we now introduce the 4 by 4 matrices 

the relations (^) and (|30|)-(^) can be reformulated as 

[Q, H] = [Q\H] = , {Q, gt} = QQ\ + QtQ = ^H- ei){H - 82) (34) 
while 

Relations (|3^-(p5|) are the ones of a quadratic deformation of the superal- 
gebra sqm{2) subtending the usual supersymmetric quantum mechanics [^. 
This quadratic superalgebra cannot be seen directly from the Dirac equa- 
tion and therefore we associate it with a hidden supersymmetry. Let us also 
and finally notice that a superalgebra similar to the one of (|3^)-(|35|) can 
also be found in the non-relativistic context when considering second order 



Darboux transformations [10|. 



6 



4 Examples 



Let us now turn to some examples and see how our method provides us new 
exactly solvable Dirac potentials from known ones. 

4.1 The free particle case 

We consider here the potential 

vo{x) = mai , X G M. (36) 

Note that it corresponds to an unusual -but convenient- realization (the 
usual one being associated with vq{x) = ma^) of the Clifford algebra sub- 
tending the one-dimensional Dirac equation. As stated i n (pO| ), it is necessary 
to take account of two eigenspinors corresponding to (|36D. Let ui and U2 
defined by 



Ul 



( cE 
' ch(A;a;) H — — sh(A;x) 

cE 

c\\(kx + 2a) + —sh{kx + 2a) 
k 



,U2 



-ch.{kx) 



(37) 



ch.{kx + 2a) 

be such eigenfunctions (of respective eigenvalues ei = E and £2 = ~E) with 

(38) 



k 



E\ 



.2a 



m — k 
m + k' 



constant. 



The unique constraint to take care of in order to apply our method is to 
have a nonvanishing determinant: detu 7^ 0. Here it is precisely given by 



detn 



1 

E 



E^c 

m + Ech {2kx + 2a) + — — sh (2kx + 2a) 

k 



4- 



(39) 



and the parameter c is such that |c| < k/E in order to satisfy this constraint. 
The result (p8[) then gives rise to the new exactly solvable potential vi 



vi{x) 



2E'c 



-0-3 



2k' 



m 



0-1 



(40) 



A " V ^ 

whose eigenspinors can be obtained from the application of L defined in 



Eq. (24) to the eigenspinors of the free Dirac Hamiltonian. Notice that the 
potential vi (x) given in Eq. ( pO| ) reduces to the well known one-soliton scalar 
potential when c = 0. 
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4.2 The generalized Coulomb case 

Before going to this example, we would like to mention that the usual radial 
equation associated to the (3+l)-dimensional Dirac equation is included in 
our developments. Indeed, the standard radial equation when coupled to 
scalar W{x) and vector V{x) potentials is 



d k 

dx X 



C73 + [M + W{x)] ai+i[E- V{x)] as | ^{x) = 0, x E M+ (41) 

where M and E are the mass and the energy of the particle while k is related 
to the total angular momentum. Eq. (pl|) can also be written as 

i<T24- + -^1 + [M + W{x)] as-lE- V{x)] ] ^P{x) = 0, (42) 
dx X J 

which coincides with hotpi^x) = Eip{x) with Hq defined in Eq. (|7|) and 

v^^{x) = ^, w?i(x) = M+V{x) + W{x), v^^ix) = -M+V{x)-W{x). (43) 



Let us now turn to our example. It corresponds to the choices of Ref. |11[| : 

V{x) = -, W{x) = ^. (44) 

X X 

We refer to this example as the generalized Coulomb one because the choice 
[a = ,B = 0) leads to the standard Coulomb interaction. 

Let ip{x) = {iIji{x),iI)2{x)) be a solution of Eq. (p|) or equivalently (p^, 
when the interactions ( ^^ are taken into account. Using standard devel- 
opments, we easily find the solutions in terms of hyper geometric confluent 
functions 



n i-Fi(l — n, 2/Li + 1; 2\nx) 



aM BEn , , , , 

-k+^ + ^] iFi(-n,2^ + l;2Anx) 



An A,^ 

n iFi(l — n, 2// + 1; 2Anx) 



(45) 



, , aM BEn , , , , 

+ [-k+^ + ^] iFi{-n, 2fi + 1; 2Anx) 

An A, 



(46) 



where the parameters An and fi are constrained by 

\l = M^-El (47) 



^2 = ^2^^2_^2 (43^ 
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while the number n is defined by 



n 



+ \- A* 



K 



This relation can be solved for the energies En as 

-a(5 ± (n + ix)-sja^ + {n + fif - /?2 



Er, 



[a2 + (n + ^)2] 



M, 



(49) 



(50) 



the plus or minus sign, as well as the values taken by n, having possibly to 
be chosen in order to ensure the square-integrability of ^'1(2^) ^-nd ^'2(2;). 
The most straightforward way to apply our method is to choose 



Ul 



ipiix) 



,U2 



n=0 



n=l 



In order to avoid heavy notations, we rewrite these choices as 



Ul 



'%1-C2X) 



with Ao and Ai defined through Eqs. (| 

^ — k 



and (pO), while 



ci 

C3 



a — /? ' 

Ai {El + M){^i-k) 

1 + 2^ (a-/3)(l + 2/i) ' 

Ai (M-gi)(a-/?) 

1 + 2^ {^l-k){l + 2^l) ■ 



(51) 

(52) 

(53) 
(54) 
(55) 



Applying finally the result (|28|), we obtain a (new) exactly solvable po- 
tential of the type 



vi{x) 



a 

- + 

X 



-M + 



+ 



k 



C2 - C3 

{ex - 82, 



C2 - C3 

X 



0-3 



+ 

X C2 - C3 



1 \ 1 



Cl 



+ 



Cl J X 



C2 
Cl 



C1C3 



ai. (56) 



In other words, we obtain a shape-invariant potential with respect to vq{x). 
A particular example corresponding to the choices 

a = 1, /3 = -1, /i = 1, A; = 1 (57) 

can be useful to illustrate the results here. Indeed we have 

Ao = 0, Ai = ^M, Cl = 0, C2 = ^M, C1C3 = ^M. (58) 
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The corresponding energies are 



El 



5 



The resulting potential is given by Eq. (56), i.e., 



1 /3M 1 

X \ b X 



f^3+ - 

X 



4M 



(59) 



(60) 



One can then determine the operator L, as defined in Eq. (|24D, connect- 
ing the eigenfunctions related to ^0(2;) = (1/x) + (M — 1/x)(T3 + (1/x)(Ti and 
to vi (x) given in Eq. (^0|) , respectively. It is given by 



L 




2M 2 

d ^4M^ 2 
d-x^ — --x 



(61) 



Due to its definition, it is clear that Lui = Lu2 = 0, with ui and U2 
of Eq. (|5^) with the values (|58|). The definition of L also implies that, 
whenever applied to any of the functions ip{x), it will give the eigenfunctions 
corresponding to vi{x) as expressed in Eq. (|60|). For instance, for n = 2, we 
have 



L 



■25^"P 



6 

125 



3M ,„ o r) 

exp 5--^ M^x^ 



50x - SOMx^ + SM^x^ 
3(-10Mx2 + 3M2x^) 

-10 + 3Mx 
5 + 3Mx 



(62) 



and one can directly check that this is a solution of the final equation 

4 

hi^jjix) = ETp(x) with E = — M. The other values of n (= 3,4,...) ev- 

5 

idently lead to similar results. 

The last information we mention here is the possibility of obtaining new 
exactly solvable potentials and not only shape-invariant ones. This situation 
arises for example when we choose 



Ul 



,U2 



n=l 



ip2{x) 



(63) 



n=2 



With the set of parameters (pT]), we obtain 

1 



vi{x) 



50x - 15Mx2 + 12M2x3 
100 + 90Mx - 
100 - 115Mx 



OOM^x^ 



100 - 115Mx 



27M\2 + 12M^x^ 



27M'^x^ + 12M^x^ 



-120Mx + 84M2x^ 



(64) 
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which has a different shape with respect to vq{x) and to the vi{x) given in 
Eq. (|60|). Evidently, one can determine the eigenfunctions related to this 
potential vi{x) defined in Eq. ( |64| ) through the application of the corre- 
sponding Darboux operator L on the eigenfunctions ip{x) of Hq. One can 
also proceed in a similar way with different values of n and obtain families 
of new exactly solvable potentials vi (x) whose eigenfunctions will be known 
through the application of the ad-hoc Darboux operator on the solutions of 
the generalized Coulomb problem. 
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